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1. Introduction 

The Hubbard model ^ is the 'standard model' for correlated electron systems 
It contains the essential ingredients for the description of itinerant electrons in a 
solid whose kinetic energy competes with their mutual Coulomb interaction. Since it 
contains very few parameters it provides an ideal development area and testing ground 
for (approximate) analytical and numerical techniques. 

Despite its simplistic structure, the Hubbard model is very difficult to solve 
exactly. An exact solution via Bethe Ansatz exists in one spatial dimension, see Ref. [3] 
for a comprehensive account. Unfortunately, an exact solution is not possible in the 
limit of infinite dimension U, and approximate numerical and analytical methods 
must be employed to solve the self-consistency equations of the resulting Dynamical 
Mean-Field Theory [3 E) • Such approximate investigations must be supported by 
perturbatively controlled analytical calculations in the limits of weak and strong 
coupling. In this way, the quality of various approximation schemes was assessed 
in Refs. |71|H1EI. 

In Ref. [7j , Feynman-Dyson perturbation theory was used to calculate the single- 
particle self-energy up to fourth order in the Hubbard interaction U. The ground- 
state energy Eq{U) and the momentum distribution n„{e{k)-,U) were obtained as 
moments of the spectral function which cannot be expressed in terms of a Taylor 
series in U . Therefore, Feynman-Dyson perturbation theory does not provide Eq(U) 
and na-{e{k);U) in form of a Taylor series in U. The coefficient of the fourth-order 
term of the ground-state energy as determined in Ref. [7] does not agree with results 
from recent Quantum Monte-Carlo investigations 

In this work we employ the Brueckner-Goldstone perturbation expansion to the 
Hubbard model in infinite dimensions. This (stationary) perturbation theory has the 
advantage that it provides Eq{U) and na{e{k); U) as a Taylor series in U. As a result, 
the Brueckner-Goldstone approach requires less numerical effort to calculate Eo{U) 
and n„{e{k); U) than the Feynman-Dyson perturbation theory and, thus, gives more 
accurate results for the same amount of computational effort. An exception is the 
quasi-particle weight Z{U) which is algebraically related to the slope of the real part 
of the self-energy at the Fermi energy and gives the size of the discontinuity of the 
momentum distribution at the Fermi energy. Therefore, we can directly compare 
Z(U) from Feynman-Dyson perturbation theory in Ref. 7 with the corresponding 
expression from Brueckner-Goldstone expansion in this work. 

Our work is organized as follows. In Sect. El we introduce the Hubbard model 
and our physical quantities of interest. In Sect. |31 we use Brueckner-Goldstone 
perturbation theory to calculate the ground-state energy of the Hubbard model in 
infinite dimensions up to and including fourth order in the Hubbard interaction. We 
find that the fourth-order coefficient quoted in Ref. |Zj is indeed incorrect. The correct 
value agrees with recent Quantum Monte-Carlo data ^Hl- In Sect. 01 we calculate 
the momentum distribution starting from the Brueckner-Goldstone expressions for 
the ground-state energy. In infinite dimensions, na-{e{k);U) can be expressed as a 
functional derivative of the ground-state energy with respect to the bare dispersion 
relation e(fc). Our results for the momentum distribution and the quasi-particle weight 
agree very well with those obtained previously from Feynman-Dyson perturbation 
theory [J]. We conclude in Sect.O 
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2. Physical quantities 
2.1. Hamilton Operator 

We investigate spin-1/2 electrons which move on a c?-dimensional lattice with L lattice 
sites. The corresponding operator for the kinetic energy reads 

f = ^ t(/ - m)c+,c„,, , (1) 

where c^^, c/^o- are creation and annihilation operators for electrons with spin a =t, J, 
on site I. The thermodynamical limit, L oo, and the limit of infinite dimensions, 
d —>■ oo, are implicitly understood henceforth. 

The kinetic energy is diagonal in momentum space 

T'-E^W^tA. (2) 

with the dispersion relation 

e{k)=ytir)e-"^^ . (3) 



Note that all momenta are taken from the first Brillouin zone. The ground state of 
the kinetic energy is the Fermi-sea |FS), 

n c+Jvacuum), (4) 

k,a-,e{k)<EF 

where Ep is the Fermi energy. 

Later we shall work with a semi-circular density of states for non-interacting 
electrons, 



^o(e) . i E - - (I)' e {iW/2f - e-)) , (5) 

k V ^ ^ 

where W = 4t is the bandwidth and 8(a;) denotes the Heaviside step function. This 
density of states is realized for a Bethe lattice with infinite connectivity In the 
following, we set t = 1 as our energy unit. 

The electron-electron interaction is taken to be purely on-site with strength U, 

H, = Ub = uY, [ni^^ - [ni^^ - ^ , (6) 

where ni^a — ^ta'^i a local density operator at site I for spin a. The Hubbard 

Hamiltonian then becomes 

H=.f + Hi=f + UD. (7) 

The Hamiltonian is invariant under SU(2) spin transformations so that the total spin is 
a good quantum number. We demand that the model exhibits particle-hole symmetry 
on bipartite lattices, i.e., H is assumed to be invariant under the transformation 

PH: c+,^{~iyci^a ; Q,. ^ (-l)'c+^ , (8) 

where (—1)' = 1 on the A-sites and (—1)' = —1 on the i?-sites. Therefore, the model 
is altogether S0(4)-symmetric jH]. This implies that there exists a nesting vector Q 
such that exp(iQZ) = (—1)' and, therefore, e{k + Q) = —e{k). We consider exclusively 
a half-filled band where the number of electrons N = 2Nr, equals the number of lattice 
sites L so that we may set the Fermi energy to zero, E-p = Q. 
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2.2. Ground-state energy, momentum distribution and quasi-particle weight 

Let |0) be the normalized exact ground state of H. As a consequence of the particle- 
hole symmetry 0, the ground-state energy of the Hubbard Hamiltonian {Tj) at half 
band- filling is an even function of U 

EoiU) = {0\H\0) = Eo{-U) . (9) 

Therefore, the Taylor expansion of Eq{U) contains even powers in U only, 

Eo{U) = + + + . . . . (10) 

For a semi-circular density of states, the average kinetic energy of the half-filled Fermi 
sea is readily calculated as 

p(o) /.O 



_ deDoie)e = ~-— . (11) 

We shall calculate the coefficients to second and fourth order with the help of 
Brueckner-Goldstone perturbation theory in section |31 
The momentum distribution is defined by 

nUU) = {0\JT.^'^'~"'^'K^^r,^JO) ■ (12) 

l,m 

In the limit of infinite dimensions, the momentum distribution depends on the 
momentum k only through the dispersion relation e(fc) 2 . Therefore, we write 
nk,a{U) = na{e{k) \ U) = nc{e; U). The particle-hole transformation ||SJ) ensures that 

na{e\ U) ^ na{e; -U) ^ 1 - n,(-e; U) . (13) 

Therefore, the Taylor expansion of n^(e] U) contains even powers in U only, 

n,{e;U)^n(°\e) + U^n^^\e) + U^n(^He) + ... , (14) 

and we may restrict ourselves to e > in the calculation of the momentum distribution. 
The leading term in H14|l is the momentum distribution for free fermions, 

6(4 -.^)e(-.).{ J t'lilil- <i^) 

We shall calculate the coefficients to second and fourth order with the help of 
Brueckner-Goldstone perturbation theory in sectional 

The discontinuity of the momentum distribution at the Fermi energy Ep = is 
of particular interest because it gives the quasi-particle weight Z(U). By definition, 

ZiU) = n^ie = 0-) - n^ie = 0+) = 1 - 2n,(e = 0+) , (16) 

and the Taylor expansion reads 

Z([/) = 1 + Z(2)^+ZW^... . (17) 

In infinite dimensions, this quantity is related to the slope of the self-energy at the 
Fermi energy via 

-1 



Z{U) 



(18) 



The self-energy E(w; U) was calculated in [7] to fourth order in the Hubbard 
interaction. Therefore, we can directly compare our results for Z{U) from Brueckner- 
Goldstone perturbation theory with those from Feynman-Dyson perturbation theory. 



Brueckner-Goldstone perturbation theory for the Hubbard model 



5 



3. Ground-state energy 

According to Brueckner-Goldstone perturbation theory |12l the ground-state 
energy can be expressed as a power series in the interaction strength, 



oo 

Eo - = {FS\Hi J2 ((4°^ - T)-'H,y |FS}e . (19) 
j=o 

Here, the subscript 'c' indicates that only connected diagrams must be included in a 
diagrammatic formulation of the theory. Note that no diagrams with Hartree bubbles 
appear because they have been subtracted explicitly in the definition of Hi. Moreover, 
due to particle- hole symmetry, only odd j contribute to the series in H19(l . 




q 

Figure 1. Second-order Goldstone diagram. 

3.1. Second order 

The ground-state energy to second order is given by 

£;^2) ^ {FS\D{E^°'^ - f)-^b\FS)c , (20) 

to which only a single diagram contributes, see Fig. ^ In this figure, a line with a 
down-going arrow indicates a particle line, a line with an up-going arrow represents a 
hole line, and a wiggiy line indicates the transferred momentum q. Note that q ^ 
applies because we can write the interaction in momentum space as 

k.p q^O 

The diagram is readily calculated from the Goldstone diagram rules, 



L 



^^[(l-n(°)(e(fc + g))]n(°)(6(fc)) 

k,p q^O 



(22) 



e(fc + q) - e(fc) -|- e{p - q) - e(p) 

The diagram was evaluated in one, two and three dimensions in Ref. 

To make progress in the limit of infinite dimensions @], we introduce the joint 
density of states via 

1 



^,(ei, £2) - ^ 2^ S{ei - e(fc + g))(5(e2 - e(fc)) , (23) 
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and write 

L J £1 - £2 + £3 - £4 

i^D,(ei,e2)i?-,(£3,£4). (24) 

We note that the joint density of states factorizes in infinite dimensions 

^(;(£ii £2) = ^o(£i)^o(£2) for 'almost aU' g. (25) 

Therefore, the contribution of the f electrons and the ]. electrons separate. This 
amounts to the statement that the second-order diagram is a 'local diagram', i.e., the 
momentum conservation at the vertices can be ignored. With the help of the Feynman 
trick 

i = / dAe-^"^ for X > (26) 
X Jo 

we write 

p{2) poo 

-jr = -J^ dA.T(A).x(A), (27) 

where 

s,(A) -Fp(A)^^h(A) (28) 
with the particle and hole contributions 

^p(^) = Z E ^"^'^"^ = I deDo{e)e-^^ , (29) 



p,e(p)>0 



^h(A) = 7 E ^ f dei?o(£)e'^ . (30) 

pAp)<o -^-"^ 

For the calculation of the ground-state energy we may set Fp{X) = Fh{X) = F{X) and 
we arrive at the final result IH 



■ dA[F(A)]^ (31) 

(2.) 

For the semi-circular density of states we find the numerical value /L = 
-0.020866148383. 



3.2. Fourth order 

The calculation of the ground-state energy to fourth order requires the evaluation of 
E^"^^ = (FS|i)(i;^"^ -f)~iL'(4°^ -f)-ii?(i;^"^ -f)-iL>|FS)c .(32) 

The Goldstone diagrams now contain four vertices and the nine diagram parts for 
one spin species are shown in Fig. |21 together with their sign which results from the 
fermionic commutation relations after the application of Wick's theorem. Note that 
the diagrams within the first, second, and third line are particle-hole symmetric to 
each other, i.e., (a) ^ {g), (b) ^ (/), and (e) ^ {h). The diagrams (c), {d) and (i) 
map onto each other under a particle-hole transformation. 
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Figure 2. Diagram parts for one spin species for the fourth-order Goldstone 
diagrams. 



In order to generate all possible diagrams, we draw all 81 possible pairs of 
diagrams from (a)| x (a)|, (a)| x to x and connect all vertices 

with horizontal interaction lines. As an example, the diagram (a)| x (z)^ is shown in 
Fig. El 

Not all of the 81 diagrams are connected, namely, the diagrams (c)| x (c)^, 
(c?)l X (d)i, and x are disconnected and drop out. Moreover, the diagrams 
{a)a X (c)_cr, (6)cr X (c)-^, {g)cr X (c)_cr, and (/)^ x (c)_cr as well as {b)„ x (d)^^ 
and (/)ct X ((i)_cr contain a hole-line and a particle-line with the same momentum and 
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therefore vanish due to the factor ria\t){\ — rLa\e)). Nevertheless, it is quite tedious 
to work out the remaining 66 diagrams even when the spin-flip symmetry is employed. 

We proceed differently. If all 66 diagrams were local diagrams, we could ignore 
the momentum conservation at the vertices as in the case of the second-order diagram 
in Sect. 18.11 The contribution of all non- vanishing connected diagrams would be 



(4) /.oo /.oo 







(-1)3 / dAi / dA2 / dA 



l2 



[a + h + c + d + e + f + g + h + i] 

- 2 [a + 5 + .g + /] c - 2 [6 + /] d , (33) 
where all functions in the integrand are functions of Ai, A2, A3. Explicitly, we have 

a(Ai,A2,A3) =Fp(Ai+A2 + A3)Fh(Ai)i^h(A2)i^h(A3) , (34) 
5(Ai, A2, A3) = Fh(Ai + A2 + A3)Fp(Ai)f^p(A2)Fp(A3) , (35) 
fo(Ai,A2,A3) = -i^p(Ai)Fh(A3)Fp(A2-f A3)Fh(Ai+A2) , (36) 
/(Ai,A2,A3)= -i^h(Ai)Fp(A3)Fh(A2 + A3)Fp(Ai+A2) , (37) 
e(Ai,A2,A3) = -i^p(Ai+A2 + A3)Fh(Ai+A2)Fp(A2)Fh(A2-KA3) , 

(38) 

/i(Ai,A2,A3) = -Fh(Ai+A2 + A3)J^p(Ai+A2)i^h(A2)^^p(A2 + A3) , 

(39) 

c(Ai,A2,A3) = Fh(Ai)Fp(Ai)Fh(A3)Fp(A3) , (40) 
d(Ai, A2, A3) = Fp(Ai + A2)Fh(Ai + A2)f^p(A2 + A3)Fh(A2 + A3) , 

(41) 

«(Ai, A2, A3) = Fp(Ai + A2 + A3)Fh(Ai + A2 + A3)Fp(A2)Fh(A2) . 

(42) 

Eq. (|33|) gives the fourth-order ground-state energy for the Wolff model 17 where the 
Hubbard interaction acts only on one site. 

However, even in infinite dimensions not all diagrams are local diagrams. First, 
there are diagrams with a 'self-energy insertion', namely (a)o- x (i)-o- as shown in 
Fig. 121 and its particle-hole transformed counterpart {g)cr x As indicated in the 

figure, two hole lines have the same momentum. Therefore, the contribution of these 
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diagrams is 

i^OO pOO pOC 

A^2{-lf dAi / dA2 / dA3[a2+52]« (43) 
Jo Jo Jo 

with 

a2(Ai,A2,A3) =Fp(Ai+A2 + A3)^^h(Ai+A3)Fh(A2) , (44) 
52(Ai, A2, A3) = Fh(Ai + A2 + A3)^p(Ai + A3)Fp(A2) . (45) 
The contribution A replaces the terms 2(a + g)i in 




Figure 4. 'Pauli-forbidden' fourth-order Goldstono diagram {d)a x {h)-„. The 
two hole-Hnes in the right part of the figure carry the same momentum. 

Second, different contributions result from the 'Pauli-forbidden diagrams'. The 
first class of Pauli-forbidden diagrams is generated by [d)^ x {h)^„, as shown in Fig.01 
and its particle-hole transformed counterpart {d)a x {e)-^. As shown in the figure, 
there appear to be two holes with the same momentum in an intermediate state. 
Recall that Goldstone diagrams which appear to violate the exclusion principle must 
be retained |12) . The corresponding contribution is 

/•OO /"OO POO 

B = 2{-\f dAi / dA2 / dA3[/i2 + e2]d (46) 
Jo Jo Jo 

with 

/i2(Ai,A2,A3) = -Fh(Ai + 2A2 + A3)Fp(A2 + A3)J^p(Ai + A2) , (47) 
e2(Ai,A2,A3) = -Fp(Ai + 2A2 + A3)Fh(A2 + A3)Fh(Ai + A2) . (48) 
The contribution B replaces the terms 2{h + e)d in (|33|l . 




Figure 5. 'PauU-forbidden' fourth-order Goldstone diagram (e)o- X The 
two hole lines in the left part of the figure carry the same momentum. 

The second class of Pauli-forbidden diagrams results from {e)^ x {^)^a^ as shown 
in Fig. and its particle-hole transformed counterpart {h)a x («)-cr, where again 
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two holes have the same momentum in an intermediate state. The corresponding 
contribution is 

/•oo /"OO (*OQ 

C = 2{-\f dAi / dA2 / dA3[/i4 + e4]i (49) 
Jo Jo Jo 

with 

/i4(Ai,A2,A3) = -Fh(Ai + A2 + A3)Fp(Ai+2A2 + A3)Fh(A2) , (50) 
e4(Ai,A2,A3) = -Fp(Ai + A2 + A3)Fh(Ai+2A2 + A3)Fp(A2) . (51) 

The contribution C replaces the terms 2{h + e)i in (|33|l . 

Altogether, the ground-state energy to fourth order becomes 



^(4) ^(4) 



L L 



A + B + C (52) 



POO poo poo 

/ dAi / dA2 / d\3[2{a + g + h + e)i + 2{h + e)d] 
Jo Jo Jo 



The single three-fold integral can be carried out on a PC after a suitable representation 
of -F(A) in (|29|1 has been found, e.g., in terms of a power series for small and large 
arguments. 

For the semi-circular density of states, we find the fairly small value e'i^^ /L = 
7.2475 • 10~^. Thus, the ground-state energy of the Hubbard model on a Bethe lattice 
with infinite connectivity reads 

' - 0.02086614838 — + 0.0000072475^ + O ' 



L 37r ■ t 

2W C/4 / 

= 0.08346459— + 0.00046384— + O [ —r 

' (53) 

up to fourth order in the interaction {W = At is the bandwidth). Similar to the 
situation in the symmetric single-impurity Anderson model |15| . the coefficient to 
fourth order is very small. For this reason, it is very difficult to extract it from 
the moments of the single-particle Green function. In fact, in Ref. 7 , the single- 
particle self-energy as a function of frequency was not calculated accurately enough 
to give the correct sign of the fourth-order coefficient. The correct value given 
in (|53() agrees very well with latest data from Quantum Montc-Carlo calculations |10|. 
^(4),QMC/^ ^ (5 ± 2) . 10-4. 

4. Momentum distribution and quasi-particle weight 

As shown in Sect. |21 in infinite dimensions the ground-state energy can be cast into 
the form 

Eo=J2<k)n{e{k)) + {0\Hi\0) , (54) 

k 

where n{e) = 2na{e{k)) is the spin-summed momentum distribution. This equation 
shows that we can obtain na{e{k)) as a functional derivative of the ground-state energy, 
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p;-2<£(p)<0 
which, together with 155|) . leads to (|15|l . 
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(56) 



Second order 

In the limit of infinite dimensions, the ground-state energy can be written as a function 
of Fp^h(A), see eq. (|^ . For e(fc) > 0, we need 



5e{k) 



^^^^^ -0 , 4(A;e)=t(A;e)F(A) , (57) 

t(A;e) = Ae~^' . (58) 

When we apply (|55|l to the expression H27() for the second-order ground-state energy 
we find for the momentum distribution to second order 

/•OC /'OC 

42)(e)=/ dAs,(A)4(A;e)= / dA Ae-^^[F(A)]3 , (59) 
Jo Jo 

in agreement with a direct calculation from Rayleigh-Schrodinger perturbation theory. 

We show the second-order contribution in Fig. 

From the value of n), '{e) at e = we deduce the second-order coefficient of the 

quasi-particle weight, Z^^) = -0.0817484. 



w 0.025 



_P 0.02 




Figure 6. Second-order contribution to the momentum distribution for e > 0. 
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4-2. Fourth order 

We can readily express the momentum distribution to fourth order in terms of the 
functions a to 02, 52, ^2, 62, and 64, and their derivatives, 

/•OO POO POO 

= / dAi / dA2 / dAg 

Jo Jo Jo 

[a + . . . + i] [a + . . . + i'] — cc — dd' — ii' 

-{a + b + g + f)c' - (a' + b' + g' + f)c 

- [b' + f + e' + h')d -{b + f + e + h)d' 

- {a' + g' + e' + h')i - (a + .g + e + h)i' 

+ {0,2 + 92 + + ei)i' + (02 + 52 + ^4 + 64)* 
+ (/i2 + e2)d' + (/i2 + e2)rf- (60) 
With the help of the abbreviation (|58|) we can express the derivatives in the form 



a'(Ai, A2, A3; 


^) 


- i(Ai + A2 + A3; e)F(Ai)F(A2)F(A3) , 


(61) 


A2, A3; 




= F(Ai + A2 + A3) [i(Ai; e)F(A2)F(A3) 








+ t(A2; e)F(Ai)F(A3) + i(A3; e)F(Ai)i^(A2)] , 


(62) 


b'{Xi, A2, A3; 




= -F(A3)F(Ai+A2)[<(Ai;e)F(A2 + A3) 








+ t(A2 + A3;e)F(Ai)] , 


(63) 


/'(Ai, A2, A3: 




= -F(Ai)F(A2 + A3)[t(A3;e)F(Ai + A2) 








+ <(Ai+A2;e)F(A3)] , 


(64) 


e'(Ai, A2, A3; 


e) 


= - F(Ai + A2)F(A2 + A3) [i(Ai + A2 + A3; e)i^(A2) 






+ t(A2;e)F(Ai + A2 + A3)] , 


(65) 


h'{Xi, A2, A3; 




= - F(Ai + A2 + A3)i^(A2) [t{\i + X2;e)F{X2 + 


A3) 






+ <(A2 + A3;e)F(Ai+A2)] , 


(66) 


c'(Ai, A2, A3; 


^) 


= F(Ai)F(A3) [t(Ai; e)F(A3) + t(A3; e)F(Ai)] , 


(67) 


^''(Ai, A2, A3; 




= F(Ai + A2)F(A2 + A3) [t(Ai + A2; e)F(A2 + A3) 






+ t(A2 + A3;e)F(Ai+A2)] , 


(68) 


i'(Ai, A2, A3; 




= F(Ai + A2 + A3)F(A2) [t(A2; e)F(Ai + A2 + A3) 






+ <(Ai +A2 + A3;e)F(A2)] , 


(69) 


12 (Ai, A2, A3 


;^) 


= i(Ai + A2 + A3; e)F(A2)F(Ai + A3) , 


(70) 


52 (Ai, A2, A3 




= F(Ai + A2 + A3) [t{Xi + A3; e)F(A2) 








+ <(A2;e)F(Ai + A3)] , 


(71) 


^2(Al, A2, A3 


;^) 


= - F(Ai + 2A2 + A3) [t(A2 + A3; e)F(Ai + A2) 








+ t(Ai+A2;e)F(A2 + A3)] , 


(72) 


e2(Ai, A2, A3 


;e) 


= - F(A2 + A3)t(Ai + 2A2 + A3; e)F(Ai + A2) , 


(73) 


^4(Ai, A2, A3 


;^) 


= - F(Ai + A2 + A3)i(Ai + 2A2 + A3; e)F(A2) , 


(74) 


e'^{\i, A2, A3 




= - F(Ai + 2A2 + A3) [t(Ai + A2 + A3; e)F(A2) 








+ t(A2;e)F(Ai + A2 + A3)] . 


(75) 



We note in passing that after some lengthy calculations we obtained the same 
expressions from Rayleigh-Schrodinger perturbation theory |16|. Thus, stationary 



Brueckner-Goldstone perturbation theory for the Hubbard model 



13 



perturbation theory is found to be applicable for the calculation of the momentum 
distribution in the thermodynamic limit. 

In the actual numerical evaluation of the momentum distribution to fourth order, 
it is recommendable to separate the contributions from the local diagrams from those 
of the self-energy diagram, Fig. and the two Pauli-forbidden diagrams. Figs. 01 
and El which are contained in the last two lines of (jSOJ. The latter contributions can 
be written as two-dimensional integrals because two energy denominators are identical 
in these diagrams and we may use the relation 

— = dAAe"^^ for x>0. (76) 



The logarithmic singularities for e — > 0^ which arise in the three diagrams cancel each 
other, and the combined integral of all terms is convergent for all e > 0. Note that in 
infinite dimensions the metallic phase obeys the exact relation 

1 - Z{U) 



n,(e^O+) 



For the first two non-trivial orders, the logarithmically divergent slope for e 
be seen in Figs. Inland [7| 
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Figure 7. Second-order contribution to the momentum distribution for e > 0. 



In Fig. \7\we show the fourth-order contribution to the momentum distribution for 
a semicircular density of states. It is seen that the correction is small for all energies. 
Our results for the momentum distribution are in very good agreement with those 
obtained from the single-particle Green function |7]. Note, however, that the Green- 
ftmction approach provides a Taylor series for the single-particle self-energy but not 
for the momentum distribution. Therefore, the numerical evaluation of Brueckner- 
Goldstone perturbation theory is more accurate than the calculation of moments of 
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the spectral function derived from Feynman-Dyson perturbation theory. This applies 
not only for the ground-state energy but also for the momentum distribution. 

A direct comparison of the results from Brueckner-Goldstone and Feynman- 
Dyson perturbation theory is possible for the quasi-particle weight. As our fourth- 
order correction we find Z^'^^ = 0.00380158. Thus, we obtain for the Hubbard model 
on a Bethe lattice with infinite connectivity 

Z{U) = 1 - 0.0817484— + 0.00380158— + O 



1 - 1.3079744 



i2 

El 



0.973204 



t4 



El 
t^ 

El. 



(78) 



These results compare favorably with those from Ref. ^. where Z2 = 1.?>Q7[V\U^ /W^ 
and Z4 — 0.969[2]f74/Vl^4 were reported. These data also agree with recent numerical 
results from Quantum Monte-Carlo calculations |10| . 



N 




Figure 8. Quasi-particle weiglit as a function of tiie interaction strength 
from GHJ. 



We show Z{U) as a function of U /W in Fig. |S1 As expected from the size of the 
coefficients in (|78|l . the fourth-order result cannot be trusted beyond U ~ W/2. This 
agrees with the result in Ref. [7]- There, the retarded self-energy up to fourth order 
violated the constraint ImE(w; J7) < for [/ > 0.6414^. As expected, fourth-order 
perturbation theory cannot be trusted beyond intermediate coupling strengths. 

5. Conclusions 

In this work we employed the Brueckner-Goldstone perturbation theory to calculate 
the ground-state energy Eq(U), the momentum distribution n^{e{k); U) and the quasi- 
particle weight Z{U) for the Hubbard model in infinite dimensions up to and including 
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fourth order in the Hubbard interactfon U. In infinite dimensions, n^{e{k);U) can 
be obtained from the ground-state energy as functional derivative with respect to 
the bare dispersion relation e(fc). A straightforward but much more tedious direct 
calculation via stationary Rayleigh-Schrodinger perturbation theory leads to the same 
expressions. 

Apart from the fourth-order coefficient of the ground-state energy, we found a very 
good agreement of our results from Brueckner-Goldstone perturbation theory with 
those obtained earlier from Feynman-Dyson perturbation theory [Jj. The corrected 
result for the ground-state energy agrees with recent Quantum Monte-Carlo data |1(J| . 
In Ref. [7], the ground-state energy and the momentum distribution were calculated 
from moments of the spectral function. Apparently, the data for the single-particle 
self-energy in Ref. |2] were not accurate enough to determine reliably the fairly small 
fourth-order coefficient in the Taylor series of the ground-state energy. 

Fourth-order perturbation theory becomes inapplicable beyond U « W/2. It is 
prohibitive to calculate the next order because the number of local diagrams alone is 
of the order of (200)^. Moreover, five-fold A-integrations appear in sixth order in U. 
It should also be clear that the sixth order would not considerably extend the region 
in U/W for which perturbation theory is applicable. Certainly, the Mott-Hubbard 
transition in infinite dimensions |2j cannot be addressed perturbatively in U /W. 
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